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Abstract 

A q-generalization of the product densities in stochastic point processes is 
developed. The properties of these functions are studied and a q-generalization 
of the usual coefficients is obtained. This for fixed q-number of particles 
coincides with the q-Stirling numbers of the second kind. The q-product 
densities are investigated using q-Poisson distribution and this shows that 
the stochastic point processes involving consistent q-generalization are in- 
herently correlated. A closely related function to q-product densities is a 
q-generalized Janossy function and a relation between the two is established. 
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I. Introduction 

In the investigation of the cosmic ray cascades, as apphcations of the 
theory of stochastic processes [2] , one deals with a stochastic variable repre- 
senting the number of particles distributed in a continuous infinity of states 
characterized by a parameter E (the energy for example). In such a study, 
Ramakrishnan [12] in 1950 introduced the concept of Product Densities which 
played a crucial role in studying a statistical assembly of particles distributed 
over a continuous parameter. In cosmic ray cascades, the number of particles 
of a given energy E or after traversing a thickness t is considered to be a 
stochastic variable. The energy E or the thickness t is taken to be a contin- 
uous variable. Thus, a stochastic variable taking discrete values, depends 
on a continuous parameter Eort. Such processes are called Stochastic Point 
Processes since the continuous variable E or t for specific values represents 
a point in the energy or thickness axis. The use of Product Densities greatly 
simplified the treatment of the stochastic variable N. Detailed account of the 
applications of product densities in various stochastic processes ranging from 
cosmic ray cascades to biology can be found in the review of Ramakrishnan 
[13], and the books by Srinivasan [14], Bartlett [1] and Bharucha-Reid [3]. 

The moments of the distribution of the number of particles (stochastic 
variable) over specified intervals {E, E + AE) can be obtained from product 
densities. The remarkable achievement of Ramakrishnan [12] is the investi- 
gation of the higher order moments in terms of product densities of higher 
degree (> 2). This led him to introduce the CJ coefficients which can be 
identified with Stirling numbers of the second kind, for fixed N, thereby 
bringing a connection to a 'number-theoretic' result of expressing m*^ power 
of a variable in terms of lower powers using falling factorials. The sum of 
CJ (s < r) over all posible values of s is the Bell's number B^^\ 

It is the purpose of this paper to consider a q-gencralization of the prod- 
uct densities by addressing a stochastic point process for q-numbers [4] dis- 
tributed over a continuos variable. Such q-numbers will be called q-stochastic 
variable. This is carried out by studying the statistical properties of a q- 
stochastic variable (q-number) with a related q-extension of the product den- 
sities. The continuous variable E oi t upon which the q-stochastic variable 
depends, is taken to be an ordinary variable. This consideration of studying 
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the statistical properties of q-stochastic variable distributed over continu- 
ous variable E, can be used to describe light wave as a stochastic variable, 
namely, the possible field realizations are stochastic. As the fields are dis- 
tributed over the space, the q-stochastic variable approach which sans the 
integer values, is apt to describe laser fight as well as thermal light [8]. 

The precise definition of q-product densities will be given in section. III. 
During the course of the study, we introduce a q-generalization of the CJ 
coefficients which become q-Stirling numbers of the second kind when the 
q-number of particles distributed, is fixed. One of the results of this classical 
analysis is a derivation of q-Bell number (See section. IV). In the literature, 
q-Bell number is investigated using q-Boson coherent states by Katriel [7]. 
As q-Stirling numbers of the second kind and the q-Bell number are the 
results of a combinatorial analysis of the q-numbers, it will be appropriate 
to obtain them using classical stochastic process as done here, rather than 
the use of semi-classical approach in [7]. Further, in this present study, q- 
product densities are investigated using probability distribution function and 
a q-generalization of the Poisson distribution is illustrated. We find a natural 
role for the q-Poisson distribution. 

The concept of product densities in stochastic processes is very briefiy re- 
viewed in section. II. In section. Ill, the q-product densties are introduced and 
their properties derived. This leads to the q-analogue of the CJ coefficients 
which for a given number of particles coincides with q-Stirling numbers of 
the second kind. They are evaluated for a few representative cases. A recur- 
sion formula for q — CI coefficients is derived. q-Bell number is discussed in 
Section. IV. The use of q-Poisson distribution in the context of q-product den- 
sities is studied in section. V and section VI is devoted to q-Janossy function 
and its relation to q-product density function. The results are summarized 
in section.VII with discussion. 

II. Brief Review of Product Density 

Let N[E(t)) be a stochastic variable (number of particles) with paramet- 
ric values < E ait. We shall suppress t hereafter. Then dN{E) is the number 
of particles in the range E and E + dE. Ramakrishnan [12] assumed that the 
probability that there occurs one particle in dE is proportional to dE, while 
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that for the occurence of n particles in dE is proportional to (dE)'^. The 
average number of particles S{dN{E)) in dE is represented by a function 
fi{E) (product density of degree one) such that 

S{dN{E)) = h{E)dE. (1) 

Denoting the probability that n particles occur in dE by P{n), we have 

P(l) = h{E)dE + 0{{dEf), 

P(0) = l-h{E)dE-0{{dEf), 

P{n) = 0{{dEY) ; n > 1. (2) 

The r*'* moment is given by 

n 

= S{dN{E)}, (3) 
where the last step follows from (2). 

The expectation value of the product of the stochastic variables dN(Ei) 
and dN{E2) has been defined in [12] as 

£{dN{E^)dN{E2)} = f2{Ei,E2)dE^dE2. (4) 

f2{Ei, E2) is the product density of degree 2 and the right hand side of (4) is 
also the joint probability that a particle lies in dEi and another in dE2 when 
dEi and dE2 do not overlap, irrespective of the number of particles in other 
ranges. When the intervals overlap, a degeneracy occurs and so in general, 

/ / £[dN{Ei)dN{E2)] = / fi{E)dE 

r-Eu f-Eu 

+ / / f2{E,,E2)dE,dE2. (5) 

t/ E£ J E£ 

Product densities of higher degree are defined by 

fn{Ei,E2r--En)dEidE2---dEn - S{dN{Ei)dN{E2) ■ ■ -dNiEn)}, {^) 
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denoting the joint probability that there hes one particle in dEi, one in dE2 
and so on where the intervals dEi,dE2, ■ ■ • dEn do not overlap, irrespective 
of the number of particles in other ranges. 

The case of particles distributed continuously in the i?-space is of 
particular interest as then is a constant. Nevertheless, the distribution in 
ii^-space follows probability laws. We have in this case. 
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(whole range) 



h{E)dE = Nf^{E)dE, 

f,{E)dE = 1, (7) 



where f^{E)dE represents the probability that any one of the N particles 
selected at random lies in dE. It follows then 

f2{E,,E2)dE,dE2 = N{N -1) f?{E,)f,{E2)dE,dE2, 
fm{Ei,---Em)dE,---dEm = N{N - 1) ■ ■ ■ {N - {m - l))f^{E,) ■ ■ ■ f^{Em) 

dEidE2 • • • dEjn, 

fN{Ei, ■ ■ ■ EN)dE, ■■■dEM = N\ fl{E^)fl{E2) ■ ■ ■ f^{Ej,)dE^ ■ ■ ■ dEj,. (8) 



The r moment of the number of particles in any finite range AE = E^ — 
Ei, S{N^^) can be represented, after taking into account the degeneracies 
in the intervals {dEi}, by 

^(N^ae) = jZC: ■■■ fs{E,,---E;)dE,---dE,, (9) 

where the coefficients are functions of r and s and do not depend on the 
ftmction /. Applying (9) to the case in which the number of particles A^ is 
fixed, using (8) and integrating over the whole range, we get 

N' = icij^, (10) 

which identifies CJ with the Stirling number of the second kind. It is to 
be noted that (9) is a general expression for the r^^ moment of A^ in terms 
of integrals of product densities of order < r and the coefficients in (9) 
become Stirling numbers of the second kind only when A^ is fixed. 



5 



Considering the case when there is no correlation between particles in any 
two different energy ranges, the product density of degree r is factored as 



fr{Ei---Er) = h{E^)---h{Er), (11) 

and (9) becomes 

^{NIe) = jlCl{E{N^E)}\ (12) 

s=l 

where the integrals are taken over the finite range AS. It is to be noted that 
(12) coincides with the moments of a variable with a Poisson distribution. To 
see this, let P(n, /\E) be a distribution function of the number of particles 
11 in /\E and G{u, AE) be its generating function, i.e., 

G{u,AE) = Y^u''P{n,AE). (13) 

n 

The r*'* moment of n is 

Sin'AE) = Y.n'P{n,AE), 

n 

where the second line follows from (13). Expanding and regrouping this we 
have, 

^{n^AE) = i:b:{G^'\^^,^E)}^^^, (15) 

where G^^\u, AE) stands for the s*^ partial derivative of G{u,AE) with 
respect to u. For a Poisson distribution. 



P{n,AE) = ,-g(n..) (g(^Ai.))- 



n! 

we have 



G{u,AE) = e^u-i)£(n^E) _ (16) 
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Then (15) becomes 

^{n^AE) = j:K{SMY. (17) 

s=l 

Comparing this with (12), we find 6^! = CJ. Thus the case of no correlation 
between particles in any two different energy ranges is governed by Poisson 
distribution. Although the result 6^ = CJ has been shown for Poisson distri- 
bution, since coefficients do not depend on /g, (15) holdsgood generally 
with Ifg replaced by CJ. 

Ill.q-generalization of Product Densities 

We shall now extend the theory of product densities to a slightly different 
type of stochastic problem, namely, a q-extension of the stochastic variable 
while keeping the continuous parameter E as ordinary. Then the statistical 
properties of a q-stochastic variable n taking discrete values 

H = (18) 

depends on the continuous parameter E and so we have a q-stochastic point 
process. When a q-stochastic variable takes value [n{E)], we need to find the 
q-number in the range E and E + dE. The q- number of particles in the range 
E and E + dE is defined to be {[n{E + dE) — n{E)]) which is just [dn{E)]. 

Now we extend the considerations in the previous section, namely, the 
probability that there occurs one particle in dE is proportional to dE and 
this is consistent with our definition given above since [1] = 1. The proba- 
bility for the occurrence of [n] particles in dE is taken to be proportional to 
[dEY. Then the average number of particles in the interval dE, denoted by 
8{[dn{E)]), is represented by a function fi^{E) such that 

£{[dn{E)]) = fi'\E)dE. 

Denoting the probability that [n] particles in the interval dE by Pg{[n]), we 
postulate (q-analogue of (2)), 

n([l]) = n(l) = fi'\E)dE + OiidEn 
PM) = P,iO) = l~fi'\E)dE-0{{dEf), 
PM) = 0{{dEY)- n>l. 
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The r^^ moment in (3) gets generalized to 

n 

= S{[dniE)]}. (19) 

The last step in the above equation follows from the postulate stated above. 
Since £{[dn{E)]) is proportional to dE, we are able to maintain the result 
of [11] that all the moments of the q-stochastic variable [dn{E)] arc equal to 
the probability that the stochastic vriable assumes the value [1] = 1. 

We shall denote the q-product densities by superscript q, 
as fl'^\E), f^^\Ei, Ei),--- fi^KEi, E2,--- E^). The q-product densities are 
necessary for the following reasons. In attempting the q-extension, one con- 
siders a distribution of [A^] particles in the energy axis. 

That is, [A^] number of particles to be distributed in the intervals dEi^ dE2, ■ ■ 
with [1] = 1 particle in each energy interval. This can be done in [^^C[i] = 
[A^]!/([1]![A/" - 1]!) = [N] number of ways, where MC[i] is the q-Binomial 
coefficient, for the first interval dEi. The remaining particles are [A^] — 1 
and also this is the number of ways of putting in the second interval. So, 
for product density of degree 2, the joint probability of putting [1] particle 
each in dEi and dE2, will be proportional to [iV]([iV] — 1) with the same 
/2 function as in (8). Proceeding in this way, with the same f function as 
in (8), we encounter [N]{[N] - 1){[N] - 2) ■ ■ ■ {[N] - N + 1) factor for the 
A^*'* interval. This will not exhaust the total number of particles. Thus, in 
dealing with q-numbers to be distributed, we need to introduce q-product 
densities such that the joint probabihty of putting [1] particle each in dEi 
and dE2 will be proportional to [A'"][A'" — 1], that for putting [1] particle each 
in dEi, dE2 and dE^ proportional to [A^][A^ — 1][A^ — 2] and so on. This will 
lead to the joint probability of putting [1] particle each in dEi, ■ ■ -dE^ pro- 
portional to [A^]!, thereby exhausting the total number of [A^] particles. The 
use of these q-product densities greatly simplifies the treatment, as will be 
shown subsequently. 

As the continuous parameters Ei, - ■ ■ En arc taken as usual variables, the 
integration overdEi ■ ■ ■ dE^ will be ordinary integrals. So, (5) and (6) become 

f-Eu rEu rEu , 

/ / S{dN{E^)dNiE2)} = / fi''\E)dE 

JEf JEi JE/ 
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fl^\Er,---En)dE^---dEn = S{dN{E^)---dN{Er,)]. (20) 

In proceeding further, the case of [N] number of particles distributed 
continuously in the £'-space is of interest. Here [N] is a constant. However, 
the distribution in E'-space follows probability laws which we are to find out. 
Therefore, we can then write in this case, 

f?{E)dE = [N\f['^\E)dE, 
I f^^\E)dE = 1, (21) 

J {wholerange) 

where fi^^^{E)dE represents the probability that any one ([1] = 1) of the [N] 
particles selected at random lies in dE. Then, 

f^'\E,,E^)dE,dE^ = [N][N -l\f['^^\E,)f['^^\E^)dE,dE^, 

f^\E,,---EjdE,---dE^ = . J^^' ,j i'^\E,) ■ ■ ■ f^^\Em)dE^ ■ ■ ■ dE^, 

[iV — my. 

ff^\E,,---Ej,)dE,---dEj, = [N]\fi''^'{E,)---fi''^''{Ej,)dE,---dEj,. (22) 



The r*'* moment of the q-number of particles in any finite range AE — 
Eu — Eg, namely ^{[-/Vj^^;} is represented by 



^{{NYae} = EC:/ •••/ fi'\E,,...E,)dE,...dE,, (23) 

g—l J Ei J El 

where is the q-analogue of the coefficient in (9). From (23) and (22), 
using (21), for [N] fixed, we obtain, 

when integrated over the whole range of dE. This expresses for [NY in 
terms of 'falling q- factorials'. The coefficients can be identified with the 
q-Stirling numbers of the second kind. By writing 
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and comparing this with the identity [A^]'""'"-^ = [A'"]''[A^] and using (24), we 
obtain a recursion relation for Cl as 



cr^ = q'-'cu + m. (25) 

Numerical values for coefficients can be obtained from (24). For r — 1, 
(24) gives Cl — 1. From the definition of [n] (18), it follows, s < n, 

[„-.] = M^. (26) 

For r — 2, using (24) and (26) we find, 

Cf ^ 1 ; C', = q. (27) 

Similarly, for r = 3, we find, 

Cf = l, C| = g([l] + [2]) ,Cl = q'. (28) 

The values in (28) can also be obtained using the recursion relation (25) and 
(27) (as input) with the conditions C^ = 0, if s = or r < s. Higher order 
Cl coefficients can be obtained recursively from (25). 

IV.q-Bell Number 

In the usual case described in Section.II, the Bell number is defined as 

= ^c;, (29) 

s=l 

and a series expansion for B^^^ is given by Dobinsky's formula [11] 

B'" = ^E^- (30) 

^ n=l "'■ 

In view of the q — coefficients CJ, we can define the q-Bell number as 

= ^c;, (31) 
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which gives using (25)- (28), 



= 1, 

= [2], 

B(^) = l + g{[l] + [2]} + 5^ 

B(^) = l + ?{[lf + [2]^ + [l][2]} + g=^{[l] + [2] + [3]} + g^ (32) 

and so on. A q-Dobinsky formula [9] can be obtained from (24). Rewriting 
(24) as 

M ^ ti^'W^- ^^^^ 

and multiplying by and summing over N from 1 to oo, we find 
Setting A = 1 we have 

ts- - - w)¥3. '''' 



s=l 



It follows from (35) that -B^^) = 1. And -B^^) = _i^^~^^i|L. Writing 

[Nf = [N] + g[A^][A^ - 1], we realize B^'^^ = [2]. Similarly it can be veri- 
fied that (32) is reproduced by (35). Although the q-Bell number and the 
corresponding Dobinsky formula have been obtained by Katriel [7] using 
q-coherent states, our derivation here uses the q-product densities and con- 
situtcs an independent approach. A rigorous derivation of the formula (35) 
for the q-Bell number using q-Poisson distribution will be given at the end of 
Section. V. 

V.q-Product densities and q-Poisson distribution 

So far, we have not used any specific distribution for q-product densities. 
We first recall a result in the theory of product densities, namely, ifG{u, AE) 
is the probability generating function to P{n, AE), the probability distribution 
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function of the total nvmber of particles n in the energy range AE and if 
fk{Ei, ■ ■ ■ Ek) is the product density of degree k, then 

U-Lh(E.r-E.)m-m ^ f^^'l^, (36) 

where the integration with respect to each variable is over the finite range 
AE. This can be derived by comparing (9) with (15) and using = b^. 

Now, we wish a q-analogue of (36). In our earher steps, we have consid- 
ered q-number of particles [n]. Prom the definition of the generating func- 
tion G{u,AE) in (13), it follows that the variable u is associated with n. 
Therefore when n is replaced by the q-number [n], the variable u should 
be treated consistently such as the differentiation with respect u should be 
q-differentiation now. Then the q-analogue of (36) will be 

where the integrations over E are ordinary integrals. The above q-analogue 

means if Q{u, AE) is the probability generating function to Pq(\n], AE), the 
probability distribution function of the q-number [n] of particles in the energy 
range AE and if fjf \Ei, ■ ■ ■ Ej.) is the q-product density of degree k, then 
(37) relates the two through q-differentiation. In order to prove (37) we 
proceed as below. The probability generating function Q{u, AE) is related 
to the probability distribution function -Pg([n], AE) by 

g{u,AE) = ^«"P,([n],AE). (38) 

n 

The r*^ moment of [n] is given by 

^{Hae} = J:HPM,^E), (39) 

n 

which can be written as 

nHAE} = {(«|r)'^(«'Ai^)} , (40) 

u=l 
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since 



OqU „ 

using (38). The terms in (40) can be regrouped to give 

^{Kae} = E^:{^^^H«,AE)}^^^, (41) 

s=l 

where G^'^^u, AE) stands for the s^^ q-derivative with respect to u. 
Now using (38) in (37) and noting 

dqu'' „ [n - fcj ! 

we have 

/ ■■■/ ft\E,,...E,)dE,...dE, = £ J^P,([n],Ai?.)(42) 

J El JEk n-k=0 F - /«J ' 

The above relation holdsgood for any q-distribution. Now, we speciahze 
in q-Poisson distribution. The generating function for q-Poisson distribution 
is 

gPoisson{u,AE) = ef (N.Ai.) ^-.([n],Ai.)^ ^^3^ 

1 

where is the q-expponential function J2'^=o ^ ^-^d ef is the second q- 
exponential function with q replaced by - and we have e^^ = 1 [10]. Then 
(38) becomes 

" n=0 n=0 

SO that the q-Poisson distribution Pj^^\[n], AE) is given by 



P^iHAE) = e7^(N,A^) £!(Ml^. (44) 
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Substituting (44) in (42), 
J El JEk [n-k\\ c, [n\\ 



= £\[n],AE). (45) 



From (23), for r = 1 and using Cl = 1, 

rEu 

Thus, 

Wi-Eu 
— i Et 

Using (47) in (23), we get for q-Poisson distribution 

£{[n]\AE} = Y^Cl{pif\E)dE). (48) 

Next, from the generating function (43) of the q-Poisson distribution, it 
follows 



E{[nlAE} = rfi'\E)dE. (46) 

J Ed 



[ ■■■[ fi'\E,,...E,)dE,...dE, = ( f^\E)dE) . (47) 

J E-i J Eu J Ep 



Q^'\u,AE)l^, = {8{[nlAE) 

= { j^^ f['\E)dE)\ (49) 



so that (48) becomes 

S{[n]\AE} = jZC:{g^'\u.^E)}^^^. (50) 

s=l 

By comparing (41) with (50), 

Cl = b:. (51) 
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As CI coefficients do not depend on (51) holdsgood in general. Now 
comparing (23) and (41) and using the above result (51), (37) is proved. 

The formula (35) for q-Bell number will now be derived using the results 
of q-Poisson distribution. The general result (23) for the specific case of 
q-Poisson distribution upon using (45) becomes 

^{Hae} = j:c:{S{[nUE}y. (52) 

Using the q-Poisson distribution (44) as 

P,in,X) = eT"^, (53) 

we have 

oo 

^{[NUe} = J2[N]P,{N,X) = A, (54) 

N=Q 

and 

oo 

^mVAE} = EWP,{NA). (55) 

Ar=o 

Then, setting A = 1, (52) leads to (35). 

VI. q-Product Densities and q-Janossy Function 

Janossy [5,6] in his study of a mathematical model for nuclear cascades, 
introduced a function, which we shall call 'Janossy Function'. This func- 
tion is closely related to the 'Product Density Function' of Ramakrishnan 

who established [13] a complete correspondence between the two functions. 
The crucial difference between the two functions is contained in their very 
definitions. Namely, if there are N particles, then fh{Ei, E2, ■ ■ ■ , Eh) is a 
product density function of degree h and fh{Ei, £"2, • • • , Eh)dEidE2 • • • dEh 
is the joint probability that there is a particle in the range dE^, a particle 
in dE2, etc and a particle in dE^, irrespective of the number of particles in 
the other ranges dEh+i,dEh+2, ■ ' ' ydE^- On the other hand, the Janossy 
function is iph{Ei, E2,- - ■ , Eh) and iph{Ei, E2, - ■ ■ , Eh)dEidE2 ■ ■ ■ dEh is the 
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probability that there is a particle in dEi, one in dE2 etc and one in dEh and 
none in other energy ranges, dEh+i,dEh^2,-''jdEN. These two functions 
are related by 

fh{Ei, £'2, • • • , Eh) = J2N=h+i (N-hy. 
lEh+i ■■■Iej,'^n{Ei,E2,---, EN)dEh+idEh+2 ■ ■ ■ dE^, (56) 

a result derived in [13]. 

We have introduced q-product density functions in (22), for [A^] particles 
distributed in dEi, dE2, ■ ■ -dE^. Integrating the third relation in (22) over 
the entire range of dEh+idEfi+2 ■ ■ ■ dEjsf, and using (21), we obtain 

{Ie,^, ■■■Ie^ fUEi, EM)dEk+idEh+2 ■ ■ ■ dEr,)dE^dE2 ■ ■ ■ dE^ ^ [N]\ 
fi'^'{E,)fi'^^'{E2) ■ ■ ■ f['^\En)dE,dE2 ■ ■ ■ dEh. (57) 

For the product f['''^'^{Ei)f[''^\E2)---f['^^\Eh)dE^dE2---dEh, we use the 
second relation in (22), so that (57) yields, 

fhi.Ei, E2,- - ■ , Eh) = 
Ieu+i lEh+2 "' ^En fliiEi, E2, - ■■ EN)dEh+idEh+2 ■ ■ ■ dEN- (58) 

We now introduce q-Janossy function ipKEi, E2, ■ ■ ■ Eh) such that 
iphlEi, E2, ■ ■ ■ , Eh)dEidE2 ■ ■ ■ dEh represents the probability that there is [h\/h 
particle in dEi, [h]/h in dE2 etc and [h]/h in dEh and none in other energy 
ranges. It is to be noted here that q-Janossy function is distinctly different 
from the ordinary Janossy function and only when q — 1 they coincide. To 
relate q-Janossy function with the q-product density function, we first write 

rhiE^,E2,---,Eh) = P^{h)rhiEuE2,---,Eh), (59) 

where the factor P'^{h) ensures that there are [h] particles in the entire energy 
range. By integrating (59) over the entire range of dEidE2 ■ ■ ■ dEh and using 
(57) with N — h, we find 

^'(^) = T^/ / ■■■/ '4^l{Ei,E2,---,Eh)dE,dE2---dEh. (60) 

\n\\ J El JEo J Eh 
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Now, fh{Ei, ^2, ■ ■ ■ , Eh) with N = h can be related to /^(^i, -E2, • • • , Eh) for 
A^, by a weighted sum of the later over /i + 1 to 00, i.e., 

00 

r,{E„E„...,Eh)\[^^^^h] - E PWA'(^i,^2,---,^.), (61) 

N=h+1 

where the in the right side of (61) stands for the q-product density of 
degree h when [N] number of particles are present. Upon identifying this 
with (58) and then using (59) with h — N,we obtain 

fhiEl,E2, - ■ ■ , Eh)\[N]=ih] = T,N=h+l [N-h]l 

Ie^+i /i?h+2 "' ^En i^N^Ei, -E2, ■ ■ ■ , EN)dEh+i ■ ■ ■ dE^, (62) 
a relation connecting q-Janossy function and q-product density function. 

V.Discussion and Summary 

We have extended the theory of product densities to q-product densi- 
ties by considering the number of particles [n] (a q-number) as a stochastic 
variable. This consideration at this stage is a mathematical possibility. It 
is also conceivable to think of q-number of particles distributed either by 
scaling the energy intervals dEi,dE2, ■ ■ ■ ,dEN and restricting to a subset 
of the intervals or by addressing a situation in which the measurement of 
the number of particles has uncertainties. This stochastic variable takes on 
discrete values (for a given q) over a continuous variable E. E is treated 
as ordinary variable and so the integrals involving E are ordinary integrals. 
Such a q-point process can be consistently described by q-product densities 
which are introduced for the first time here. The r*^ moment of [n] is de- 
scribed in terms of q-generalized coefficients with falling factorials. For a 
fixed [N] , CI coefficients are identified with q-Stirfing numbers of the second 
kind. A recursion relation for CI is obtained. This leads to g-Bell number 
and 5-Dobinsky formula which is derived here using only the q-product den- 
sities and this agrees with the result obtained using q-boson coherent states. 
The notion of q-product densities is further investigated using q-probability 
distribution and its generating function. The case of q-Poisson distribution 
and its generating function is discussed in detail. 

A q-generalization of Janossy function is given and its relation with the 
q-product density function is derived. 
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In this approach, since the continuous variable E is taken to be the ordi- 
nary variable, the calculus of differentiation with respect to and integration 
over E remain unaltered. Consequently the evolution equations for the q- 
product densities will be very similar to those of the usual product densities. 
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